Maximum bounds on the surface redshift of anisotropic stars 
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It is shown that for reahstic anisotropic star models the surface redshift can not exceed the values 
3.842 or 5.211 when the tangential pressure satisfies the strong or the dominant energy condition 
respectively. Both values are higher than 2, the bound in the perfect fluid case. 
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I. INTRODUCTION 

CN ; 

, It is well known that the surface redshift s for a static perfect fluid sphere, whose density is positive and not 
increasing outwards, is not larger than = 2 The bound holds for the interior Schwarzschild solution with 
constant fluid density p and infinite central pressure. It occurs before the appearance of a horizon. 
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' During the years, different arguments have been put forth for the existence of anisotropy in star models. It may 
' be due to the presence of solid core, phase transitions, mixture of two fluids, slow rotation The idea that the 
, tangential pressure t may be different from the radial pressure p was suggested first by Lemaitre in 1933 He 
' discussed a model sustained solely by tangential pressures and with constant density. This model was generalized for 
variable density by Florides and was recently revisited in Ref. The application of anisotropic fluid models to 
' neutron stars began with the pioneering work of Bowers and Liang | and was done both analytically and numerically 
^ [0-|l5|. Some recent work may be found in Refs. [p^-pO[. 

In many papers it is stressed that arbitrarily big redshifts are obtained when t grows to infinity [5|-|^, [l0| , [l^ , p^ . 
However, in realistic models t should be finite, positive, and should satisfy the dominant energy condition (DEC) 
t < p 01 even the strong energy condition (SEC) 2t + p< p. These may be written together a.s t < ep where e = 1 for 
' DEC and e — 1/2 for SEC (if the realistic condition for positive p in the interior is accepted). 

, The bounds on s in anisotropic models were studied in Ref. [ pl| in a more general setting which incorporates soap 
■ bubbles, monopoles and wormholes and the focus was that a horizon does not form. It generalized this result in the 

Q ' perfect fluid case |^ to anisotropic fluids. Similar conclusions were reached in Ref. [ p3| . 

In this paper we elaborate further on the method used in Ref. pi]] and give concrete values for the maximum surface 

^ I redshift when the tangential pressure satisfies either DEC or SEC. 

bJO In Sec. II the field equations are given in a convenient form and the bound on t is implemented to derive an 

^ inequality for the mass-radius ratio. 
• ^ In Sec. HI this inequality is used to obtain maximum bounds on the surface redshift when SEC, DEC or a 

pseudoisotropy condition is satisfied by t. Several realistic models are discussed and it is shown that the bounds can 

H , not be saturated. 
■ - - ■ Sec. IV contains short discussion. 

II. FIELD EQUATIONS AND THE MAIN INEQUALITY 

The metric element in curvature coordinates is 

ds^ = e"dT^ - e^dr^ - [dO^ + sin^ Od^^) , (1) 
where z/, A depend only on the radial coordinate r. The Einstein equations read |^ 

np= — (x'--)+^, (2) 



* E-mail: boyko@inrne.bas.bg 



1 



«P = — (j^' + -)-\, (3) 
r V r J r'' 



^t = L-'l^"-'^ + i- + ) , (4) 



2" V' '2~^~2 ' r 
where ' means derivative with respect to r, k = SwG/d^ and we use units with G = c=l. Eq. (2) integrates to 
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2m 1^ r 2, fr\ 

1 , = 77 / pr dr. [oj 



Here m is the mass function, m > and consequently z < 1. Eqs. (3)-(4) give a linear second-order equation for 
y = e'l'^: 

2r^zy" + [r^z' - 2rz) y' + [2 (1 - + rz' - 2KAr'^] y = 0. (6) 
Here A = t — p is the anisotropy factor. Finally, Eq.(3) may be written as 

For stability reasons p must be positive and hence y' > 0. Eq. (6) has a more compact form 

.'"(^)=Dy, D=(^y + ^. (8) 

It is convenient to introduce the average density 

<P>=-,1 Pr'dr=—,. (9) 

Then Eqs. (5), (8), (9) give 

D = ^{p-<p> +2A) . (10) 

A realistic requirement is that p should be finite and positive. It must decrease monotonically or stay constant for 
stability reasons, p' < 0. Then it is easily shown that < p >'< and p << p>- Written in another way 

^ < 3. (11) 

Now, let us divide Eqs. (7) and (8): 

rfln^^V = ^-<^>+f^. (12) 

\ r J p+ < p > /3 ^ ' 

There are two cases; D < everywhere and D > somewhere. In the first case the r.h.s. of Eq. (12) is bounded from 
above by zero. Perfect fluids (A = 0) form a subcase of this case. Anisotropic fluids with A < 0, which have radially 
dominated pressure, form another subcase. Even when A > in some regions there is a third subcase with D < 0. 
As was mentioned in the introduction, a realistic t satisfies the inequality t < ep. Then A. < t < ep and a suSicient 
condition for non-positive D is 

dlnm 3 , , 

rflnr - Y+Ye' ^ ' 

When e = we return to the isotropic case and Eq. (11) although t = p instead of f < 0. Eq. (13) shows that when 

the slope of m in a logarithmic scale is not steep enough, no positive A is able to ensure I? > 0. We shall show in the 
following section that models with D < satisfy the Buchdahl bound on the redshift, while models with D > Q can 
develop in principle bigger redshifts. 
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Anisotropic models are subjected to three field equations but possess five characteristics z, y, t,p, p. Therefore, two 
of them or their combinations must be given explicitly. Eq. (13) is important for models with a given density profile. 
This profile should break at least once the criterion in order to possibly achieve redshifts bigger than 2. Another 
possibility is to satisfy directly the condition D > 0, which is more general, but can be checked without solving the 
field equations only when A is the second given function. 

When 13 > in some region, the following chain of inequalities holds 

P^<P>+2A ^ ^(t-p)^ < ^ < fe. (14) 
p+ < p> /3 p+ < p > /3 < p > 

The isotropic case is regained again when e = 0. Inserting the bound from Eq. (14) into Eq. (12) and integrating 
from r to the boundary of the fluid sphere at r = i? yields 

2/'>^W(^)"^, Mr)-'-^. (15) 

At R the interior fluid solution should be matched to the exterior Schwarzschild solution, which requires z (R) = 
1 — 2M/R. Here M = m (R) is the total mass of the fluid sphere. Eq. (9) provides the inequality 

2m 2M 2 

- > -J^r^- (16) 

After these remarks, let us integrate Eq. (15) from the centre to the boundary and take into account that y (0) > 0. 
The result reads 



^ r) -2i?3+6ey^ (1_2M^)V2- (17) 

Eq. (17) is a particular case of Eq. (66) from Rcf. In this paper polar Gaussian coordinates were used, the bound 
in Eq. (14) was assumed to hold for some positive e, not connected in general with t, and the maximum of 2m/r was 
not obliged to be on the surface of the configuration. Eq. (17) can be written as an inequality just for a = 2M/R 

"'O [1 — X) ' 

This is the main inequality to be used for finding redshift bounds. 

III. BOUNDS AND MODELS 

Eq. (18) provides maximum values for the mass-radius ratio a and the surface redshift 

s = (l-a)"'/'-l (19) 
both of which depend on e. In the perfect fluid case formally e ~ and Eq. (18) becomes 

3(l~a)'/'>l. (20) 

This gives the values found by Buchdahl, am = 8/9, Sm — 2. When anisotropy is present and t is untied from p, e > 
and there is possibility for higher redshifts. In general, the integral on the r.h.s. of Eq. (18) is expressed through 
the hypergeometric function. When Se is an integer it becomes a sum of different powers of 1 — a. Thus, when DEC 
holds for t {e — 1) we have 



(1 - a) 
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> - (21) 
- 35 ^ ' 



and the maxinmm bounds are a,n — 0.974 and = 5.211. These values were found in Ref. pl| . 

An useful relation which follows from the field equations is the Tolman-Oppenheimer-Volkoff (TOV) equation 



3 



2rz r 

It shows that for reahstic and, hence, finite at the centre p,p' we have A (0) = 0. Furthermore, the radial pressure 
is obhged to vanish at the boundary. If the more imperative SEC holds for t , p > p + 2t gives t < and 
even t (0) < ^p(O). Let us accept the stronger bound, which holds in the interior of perfect fluids, and call this a 
pseudoisotropic case, with t,p being of the same magnitude, although not equal. Then e = 1/3 and Eq. (18) becomes 

+ >1, (23) 

which yields a^n — 0.946 and s„j = 3.310. These values are still above the Buchdahl ones. 

The general value of e for SEC, however, is 1/2 and this case is the most realistic one. Now 3e is not an integer, 
but the integral in Eq. (18) is once more expressed in terms of elementary functions and gives 

4 l-a^/'a3/2 > arcsin 2a-l +T77- 24 

4\/a - 8 16 

Computer calculations show that «,„ = 0.957 and s„i — 3.842. The last number is the main result of this paper. It 
shows an almost double increase in the Buchdahl bound when anisotropy is allowed and t satisfies SEC. The other 
assumptions made were p > 0,p > 0, p' < and the inevitable t > p. Of course, D > should also be true, otherwise 
the bound collapses to the Buchdahl's one. Thus, realistic anisotropic star models can possess higher redshifts than 
the isotropic ones but it is limited and never reaches infinity. 

It is said sometimes for perfect fluids that the Buchdahl bound is optimal because there is a model which saturates 
it. Eq. (12) shows that for such a model p —< p > which is possible only for constant p. Thus we arrive at the 
Schwarzschild incompressible sphere with the central pressure as a free parameter. Only the model with p (0) = oo 
saturates the bound. A possible explanation is that y (0) > and the equality is required for saturation. This leads to 
singular metric and it is probably induced only by singular central pressure. When the pressure is required to satisfy 
SEC one has only s,„ = 1/2 We leave aside the problem that constant density leads to infinite speed of sound 

V — (dp/dp)^^^ , while a realistic speed of sound is bound by 1, the speed of light in our units. In conclusion, there is 
no realistic model, saturating the Buchdahl bound. 

What is the situation for anisotropic models? According to Eq. (12) a saturating model must satisfy 

p- < p> +2A = 6e(^p+^< p>^ . (25) 

However, this is not possible at the center because there the l.h.s. vanishes, while the r.h.s. is positive. Again the 
explanation may be connected with the requirement p{0) — oo because of y (0) = 0. Then A (0) < and Z? (0) < 
since t (0) is finite, but this contradicts the nature of the saturation {D > 0). 

One may ask for models which satisfy the limiting assumptions made during the derivation of Eq. (14). Namely, 
p = const = pq, p = and t = epo- These are three conditions while only two of the fluid's characteristics can be 
fixed beforehand. If we take the first two conditions, we arrive at the Lemaitre model [||-D. Eq. (22) then gives 

t^^Po. (26) 

Obviously t is not as required. Using the bound on t gives 

4e 

This inequality yields a„i = 0.8, Sm — 1.236 for DEC and a,„ = 0.667, Sm — 0.732 for SEC. These values are higher 
than the ones of the realistic Schwarzschild interior solution, but faraway from the absolute bounds derived above. 
Two models possessing the values when DEC holds were given recently [|l^ (models I and IV). 

Most successful is the Bondi model ||l^ which also has a constant density, but the second given function is the 
constant ratio Q = {p + 2t) / pQ. The region Q < 1 is studied, which is equivalent to SEC. The redshift increases with 
Q and numerical simulation gives Sm = 1.352 for Q 1. Close to this result stands Example 2 from Ref. |lj] with 
s — 1.2. It utilizes a nonlocal equation of state and a density profile used in Refs. Pj9|,|l^. 

Finally, the conformally flat anisotropic models are worth being mentioned The vanishing of the Weyl tensor 
implies 
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The two terms in D are of the same nature, but do not compensate each other. All characteristics of the solution 
depend only on the density profile. D is definitely positive, however, already the Buchdahl bound is reached with 
infinite central radial pressure. This is analogous to the behavior of the interior Schwarzschild solution which is the 
only conformally flat solution in the perfect fluid case Q . 



IV. DISCUSSION 



The bounds found in this paper generalize the Buchdahl bound to anisotropic star models and show that claims for 
arbitrarily large redshifts are not realistic. The bounds are absolute, i.e. numbers depending on few simple realistic 
requirements. They do not depend on the details of the mechanism generating anisotropy, equations of state, central 
density and so on. We have not explored also the consequences of the condition < -y < f , necessary for the causal 
behavior of the fluid. In the perfect fluid case a lower bound s„i = 0.854 was established heuristically when some of 
the above mentioned features are taken into account p5[ | . Probably, such a study can be performed in the anisotropic 
case too. It may concern not only the surface redshift but also the maximum masses and moments of inertia of 
neutron stars. 
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